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On Symmetry Properties of Quaternionic Analogs of Julia Sets 
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B. I. Stepanov Institute of Physics of the National Academy of Sciences of Belarus, 
F. Skariny Ave. 70, 220072 Minsk, Belarus 

By means of theory group analysis, some algebraic and geometrical properties of 
quaternion analogs of Julia sets are investigated. We argue that symmetries, intrinsic 



to quaternions, give rise to the class of identical Julia sets, which does not exist in 
complex number case. In the case of quadratic quaternionic mapping X^+i = X^ + C 
these symmetries mean, that the shape of fractal Julia set is completely defined by 
just two numbers, Co and |C|. Moreover, for given Co the vector part of the Julia 
set may be obtained by rotation of a two-dimensional Julia subset of arbitrary plane, 



comprising C, around the axis n = C/|C|. 



I. INTRODUCTION 



> 

o 
o 

Beautiful and unexpected fractal properties of Julia-Fatou and Mandelbrot sets, spawned 



by a simple iteration rule 

z k+1 = z\ + C, (1) 



with z and control parameter C being complex numbers, have been intensively studied by 
many authors [ffl, 0, [J. Iterations of arbitrary starting point zq in accordance with (II]) 
ultimately result in the confinement of all Zk,k > N max » 1 in a finite region, a basin of 



attraction, of some Zj. These Z\ are called attractors; they are completely defined by control 
parameter C. 

A key concept of fractal set research is the special case of z = oo attractor. In fact, each 
C of (|I|) classifies all points of the complex plane as belonging to either runaway subset or 
prisoners one. If z$ belongs to runaway subset, then ([]]) leads to z = oo. In the other case, 
ultimate cycles of z^ will reside in a basin of some finite attractor. 

A fascinating boundary of the basin of attraction of z = oo is called Julia set. All 
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information about shapes and topologies of these sets is encoded in C. These sets are 
mostly fractals. 
Prisoners of 

Ck+i = C\ + Ck (2) 

(i.e. of quadratic mapping ([I]) with starting point zq = 0) form the famous Mandelbrot set. 
It classifies all Julia sets as either connected or unconnected ones. The latter are known as 
Cantor dust or Fatou sets. 

There were many attempts to apply fractal sets to various fields of physics (see [|], |5j). 
It was clear from the very beginning, that many attractive features of these sets were due 
to remarkable properties of algebra of complex numbers. This suggests to generalize the 
quadratic mapping ([I]) to different algebras and to study their behavior. 

In paper || such a generalization was undertaken for double numbers, which differ from 
complex ones just in definition of the imaginary unit e, i.e. e 2 — 1. Replacement of complex 
numbers by double numbers actually means a transition from the Euclidean plane to the 
pseudoeuclidean one, the latter being of special interest for physics, in particular for rel- 
ativistic kinematics. Another example of validity of double numbers give two-dimensional 
relativistic models, which are prevalent in modern field theories. They also may be success- 
fully described in terms of double numbers 0. 

Generally speaking, all hypercomplex algebras, i.e. algebra of double and dual numbers, 
of quaternions, biquaternions and of octanions, have proved to be very convenient in nu- 
merous physical applications. It may be explained by their close relations with geometries 
of Euclidean and pseudoeuclidean spaces and with spaces of constant curvature . 

In this paper we focus on quaternion algebra. Really, the 'quaternion language' seems 
to be especially natural for physics of our four-dimensional space-time. It has been already 
successfully exploited in describing of rigid body motion ||, in searching for instanton so- 



lutions of Yang-Mills equations [TO], in the problem of supersymmetric oscillator [IT| and 
many others. 



II. QUATERNIONS 



Let us first briefly remind the fundamentals of quaternion algebra. The quaternion X 
is a set of 4 real numbers xq, x\, X2, X3 with 3 imaginary units i,j,k. The definition and 



commutation properties follows: 

X = x + ixi + jx 2 + kx 3 , (3) 

where 



(4) 



i 2 


= J 2 = 


k 2 = -l 


ij 


= -ji = 


k- 


ki 


= —ik = 


j; 


jk 


= -kj = 





Hence quaternions are non-commutative with respect to multiplication, so that 

X 1 -X 2 ^X 2 -X 1 . (5) 

In some applications a representation of quaternion (^) as a pair of complex numbers, 

X = (x + ixi) + (x 2 + ix 3 )j, (6) 

may turn to be useful. The traditional notation of arbitrary quaternion (|3|) as a combination 
of scalar xq and three-dimensional vector x parts, 

X = x + x, (7) 

may be very convenient in many problems. In terms of (0), the operation of quaternion 
conjugation is simply 

X = Xq — x. 

The product of two arbitrary elements of the quaternionic algebra (i.e. of two four- 
dimensional vectors) in the traditional notations reads: 

XX' = x x' - (x ■ x') + x x.' + x'qK + [xx'j. (8) 

Consequently, the square of any quaternion X looks like 

X 2 = XX = x 2 - |x| 2 + 2x x. (9) 

It should be noted, that operation of division by nonzero element (q ^ 0) is well defined 
for quaternions. 

The set of quaternions forms a group with respect to the operations of addition and 
multiplication. 
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III. QUATERNIONIC ANALOGS OF FRACTAL SETS 

The above-mentioned algebraic properties of quaternions permit to introduce a quater- 
nionic analog of the original complex number Julia-Fatou algorithm ([]]). Namely, 

X k+1 = X 2 k + C, (10) 

with X and control parameter C being quaternions. This iteration rule maps the four- 
dimensional Euclidean space on itself. Direct computer experiments justify, that (|10D classi- 
fies again all quaternions with respect to X = oo as belonging to either prisoner or runaway 
subsets. The border of these subsets is a quaternionic analog of Julia set; it is fully defined 
by coefficients of mapping rule, e.g. C in case of ([K]). However, in contrast to complex and 
double numbers, these fractal sets are four- dimensional. 



The first treatment of this problem was given by A. Norton in Ref. fl2fl . The author have 
discussed some aspects of quadratic mappings 

X k+1 = AX 2 k + B (11) 

with A, B and X being quaternions. It was noted, that the general case the problem was 
very difficult both for discussion and for representation of results. 

To simplify the problem, the consideration of Ref. [|12 was restricted to a special case 



of coefficients A and B belonging to the complex number subalgebra of quaternions. The 
author argued, that this choice allowed to reduce the dimensions of the resulting Julia sets. 
However, no good arguments were proposed to support this statement. 

A particular attention of Ref. |l!| had been paid to the analysis of topological properties 
of quaternionic Julia sets. It was noted, that, due to noncommutativity of quaternions, 
there actually exist three different non-equivalent generalizations of ([[]). Besides (|11|), one 
is also to discuss 

X k+1 = X k AX k + B, (12) 
Xk+i = X\A + B. 

The main goal of our research is to provide an instrument for establishing of common 
features and differences of fractal sets realized by algorithms (|T0l)-(|T^) from algebraic, the- 
oretical group point of view. 
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Let us consider transformations connected with multiplications and additions of quater- 
nions as transformations of the group of motion of four-dimensional Euclidean space || [13 
It is easy to check that (10)~(12) are invariant under the following transformations: 



X' k = QX k Q, A' = QAQ, B' = QBQ, C' = QCQ, (13) 
where quaternion Q satisfies the condition 

QQ = 1, (14) 

i.e. Q = Q-\ 

Note, that due to commutativity of complex and double numbers the analogous trans- 
formations, q = e iip , mean just trivial rotation by angle <p of the whole plane. 

Transformations (|13|) are inner automorphisms of the division ring of quaternions. They 
are isomorphic to the transformations of the group of three-dimensional rotations SO(3.R). 
Note, that Xq and |x| stay invariant under such transformations. 

Thus, we actually deal with a class of equivalent iteration rules 

X' k+1 = A'(X' k f + B', (15) 



which are bound with ( |Tl~D by relations (|I3|-|T4]) . Applying these transformations to the 
mapping (pi]), we obtain equivalent rules 

X' k+ i = {X' k f + C. (16) 

To illustrate this symmetry manifestation, let us show, that there is a freedom in ori- 
entation of vector C (remember, that C = Cq, |C'| = |C|). These transformations are 



analogous to the plane transformations of the Lorentz group, found in Ref. [|14}| . For vector 
parts of quaternions C and C the plane transformation quaternion Q and its conjugate Q 
are: 

(C + C') 2 v /z C^ 



Q - C C + ° - C + C ' °' ■ (18) 
v^C^ / (c + C') 2 J-(C + C') 2 



It may be convenient to rotate vector C so that to orient vector C along, say, vector i: 

C = C + |C|z. (19) 
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If all starting points X of ( |TTD lie in this complex plane, the following points Xj, will lie in 
this plane too. Thus, one obtains a complex number Julia set, which is a subset of the full 
quaternionic one. 

The most important example of the profit of ( |T3| ) in studies of quaternionic Julia set 
symmetries is given by transformations, that leave control parameter C intact, i.e. 

C' = QCQ = C, QQ = 1. (20) 

It is straightforward to check that the proper quaternion is (see Ref. ||): 

1 + ntan % c? 09 . . 

Q= ; =L = cos- + nsin-, 21 
" '1 + tan 2 2 2 2' 



2 

where n = C/ 1 C | . 

These transformations form the 0(2) group. They are similar to the gauge symmetries 
of field theories. One can consider the resulting Julia set as a projective space and bundle 
manifold. An invariance of fractal sets under transformation (|21| ) means that projection of 
this sets in the three dimensional subspace are axial symmetric. 

These symmetries does not change the zero quaternion component. They do exist for 
each Co, but corresponding two-dimensional Julia subsets are different even for the same 
|C|. It fact, quaternionic symmetries mean, that one needs just two numbers, Co and |C|, 
to describe all possible shapes of quaternionic Julia sets. 

Besides continuous symmetries, quadratic mappings of all hypercomplex algebras also 
possess a special obvious invariance under discrete reflections, 

X «-> -X. (22) 

This invariance is due to two- valued property of X 2 . Really, if some point Xj, belongs to a 
Julia set, then the same is true for — Xk- This symmetry allows to save the computer time 
at calculations of fractal pictures. 



It may be interesting to mention, that the special choice of parameters A and B in Ref. JT2 
as complex numbers was not obligatory. Symmetries (|T3|) imply, that it would suffice just 
to orient vector parts of both quaternions in the same direction. The resulting Julia sets 
would be obtained by rotation of arbitrary plane Julia subset around the axis oriented along 
the direction of their vector parts. The transverse cross-sections of these Julia sets consist 
of concentric circles. 
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IV. CONCLUSIONS 

Quaternion algebra provides non-trivial generalizations of usual Julia sets. Although 
these sets are much more complicated, the intrinsic quaternionic symmetries allow to simplify 
the problem. It turns out that in case of quadratic mapping (pi]) all essentially different 
quaternionic analogs of Julia sets may be enumerated by just two numbers, Cq and |C|. Due 
to 0(2) symmetry, the three-dimensional part of quaternionic Julia sets may be restored by 
rotation of some arbitrary two-dimensional (e.g. complex number) Julia subset around the 
axis n = C/|C|, that lies in this plane. 

This work was supported by the Belarusian Republic Foundation of Fundamental Re- 
search, project No. P-229/585. 
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